Waveform-Agile  Tracking  In  Heavy  Sea  Clutter 

Sandeep  P.  Sira,  Antonia  Papandreou-Suppappola,  Darryl  Morrell^,  and  Douglas  Cochran 
SenSIP  Center,  Department  of  Electrical  Engineering,  Arizona  State  University,  Tempe,  AZ  85282,  USA 
t Department  of  Engineering,  Arizona  State  University  at  the  Polytechnic  Campus,  Mesa,  AZ  85212,  USA 
Email:  {ssira,  papandreou,  morrell,  cochran}  @asu.edu 


Abstract — The  detection  and  tracking  of  small  targets  on  the 
ocean  surface  is  a  challenging  problem  due  to  low  signal-to- 
clutter  ratio  (SCR)  that  results  from  low  grazing  angles  and 
high  sea  state.  Recent  advances  in  sensing  technologies  have 
enabled  waveform-agile  schemes  that  tailor  the  sensor  waveform 
to  match  the  overall  sensing  objective.  In  this  paper,  we  propose 
a  methodology  to  dynamically  adapt  the  transmitted  waveform 
to  improve  the  tracking  performance  for  scenarios  characterized 
by  heavy  sea  clutter.  Employing  the  compound-Gaussian  model 
for  sea  clutter,  we  develop  an  algorithm  for  online  design  of  a 
phase-modulated  waveform  that  improves  the  SCR  in  a  range 
bin  of  interest.  A  particle  filter  tracker  uses  the  measurements 
obtained  by  this  waveform  to  estimate  the  target  state.  We  present 
a  simulation  study  to  demonstrate  that  our  scheme  leads  to 
improved  tracking  performance. 

L  Introduction 

Waveform  diversity  is  fast  becoming  one  of  the  most  im¬ 
portant  methods  by  which  sensing  systems  can  be  dynamically 
adapted  to  their  environment  and  task  to  achieve  performance 
gains  over  non-adaptive  systems.  With  recent  advances  in 
sensor  hardware,  this  adaptation  can  involve  the  changing 
of  the  transmitted  waveform  ‘on-the-fly’  so  as  to  obtain 
information  that  optimally  contributes  to  back-end  sensing 
algorithms  such  as  trackers  and  detectors.  The  benefits  of  such 
adaptation  include  improved  performance  and  reduced  sensor 
usage  leading  to  greater  system  efficiency.  Waveform-agile 
sensing  has  recently  received  wide  attention  with  a  number 
of  investigations  for  target  tracking  [  1  ]— [3] ,  detection  [4]-[6], 
and  classification  [7],  [8],  reported  in  the  literature. 

Most  investigations  of  waveform-agile  tracking  have  fo¬ 
cused  on  the  exploitation  of  the  delay  and  Doppler  resolution 
properties  of  transmitted  waveforms  to  estimate  the  range 
and  velocity  of  a  target  in  accordance  with  the  dynamically 
changing  uncertainty  of  the  tracker  [  1  ]— [3] ,  [9].  Here,  the 
impact  of  a  particular  waveform  on  the  measurement  errors 
was  quantified  by  the  Cramer-Rao  lower  bound  (CRLB), 
which  is  derived  from  the  curvature  of  the  waveform  ambiguity 
function  at  the  origin  in  the  delay-Doppler  plane  [10].  Since 
this  characterization  ignores  the  sidelobes  of  the  ambiguity 
function,  it  does  not  seem  appropriate  in  situations  involving 
low  signal-to-noise  ratio  (SNR)  or  low  SCR. 

Furthermore,  recent  work  on  dynamic  waveform  adaptation 
for  target  tracking  has  often  assumed  perfect  detection  [1], 
or  simplistic  clutter  models  [3],  [9],  [11].  While  this  may 
be  appropriate  when  the  SNR  is  high,  neither  assumption 
is  justified  in  scenarios  that  involve  heavy  sea  clutter.  Early 
investigations  [12]  of  the  design  of  waveforms  for  clutter 


rejection  assumed  Gaussian  models  for  sea  clutter  echoes. 
Such  models  are  now  known  to  be  inadequate  due  to  the  higher 
spatial  resolution  offered  by  modern  radars.  As  a  result,  the 
compound-Gaussian  model  for  sea  clutter  has  gained  wide 
acceptance  [13],  [14],  Although  there  has  been  extensive 
research  in  the  development  of  detection  strategies  for  targets 
in  non-Gaussian  sea  clutter  [15],  very  few  applications  of 
dynamic  waveform-adaptation  to  improve  target  tracking 
performance  in  the  presence  of  non-Gaussian  interference  have 
been  reported  to  date. 

In  this  paper,  we  approach  the  challenging  problem  of 
tracking  a  small  target  in  heavy  sea  clutter  with  the  premise 
that  improved  detection  performance  also  leads  to  improved 
tracking  performance.  This  is  due  to  the  fact  that,  for 
a  given  probability  of  detection,  a  lower  probability  of 
false  alarm  implies  less  uncertainty  in  the  origin  of  the 
measurement  which  leads  to  lower  tracking  error  [16].  Using  a 
previously  developed  procedure  for  dynamic  waveform  design 
for  detecting  targets  in  low  SCRs  [6],  we  develop  a  particle 
filter  based  tracker  that  uses  the  measurements  obtained  by  the 
adapted  waveform.  A  simulation  study  based  on  parameters 
derived  from  real  sea  clutter  measurements  demonstrates  that 
our  approach  provides  significant  reduction  in  the  tracking  root 
mean  square  error  (RMSE)  when  compared  to  a  non-adaptive 
system. 

The  paper  is  organized  as  follows.  In  Section  II,  we  describe 
the  models  for  the  target  dynamics,  clutter,  and  observations. 
Section  III  discusses  the  dynamic  waveform  design  algorithm 
and  the  estimation  of  clutter  statistics.  In  Section  IV,  we 
present  the  target  tracker  in  detail,  while  in  Section  V,  we 
provide  the  results  of  a  simulation  study. 


II.  Target,  Signal,  and  Clutter  Models 

Let  Xfc  =  [rfc  fk]T  define  the  state  of  a  target  at  time  k, 
where  and  A  are  the  range  and  radial  velocity,  respectively, 
with  respect  to  a  fixed  radar.  We  assume  a  linear  dynamics 
model  for  the  target  motion  so  that 

Xi  =  FX*_1  +  Wk)  (1) 


where  is  a  zero-mean,  white  Gaussian  noise  sequence 
with  covariance  matrix  Q.  The  constant  matrices  F  in  (1)  and 
Q  are  given  by 
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respectively,  where  AT  is  the  interval  between  two  observa¬ 
tions  and  q  is  a  constant. 


A.  Received  Signal  Model 

We  consider  a  medium  pulse  repetition  frequency  (PRF) 
radar  that  transmits  a  burst  of  2 L  pulses  in  each  dwell  in 
an  area  or  direction  of  interest  as  shown  in  Fig.  1.  Each 
dwell  consists  of  two  sub-dwells,  Sub-dwells  1  and  2,  during 
which  L  identical  pulses  of  the  waveforms  Si(t)  and  S2(t) 
are  transmitted,  respectively.  At  the  end  of  Sub-dwell  2,  a 
measurement  is  provided  to  the  target  tracker.  While  s\(t)  is 
a  fixed  waveform,  *2  ft)  is  dynamically  designed  to  improve 
the  SCR  at  a  predicted  target  location,  thus  improving  the 
tracking  performance. 

With  s(t)  denoting  either  si(t)  or  S2(t),  the  received  signal 
at  the  fth  pulse,  l  =  0, 1, . . . ,  L  —  1,  in  each  sub-dwell,  is 

g\t)  =  bls(t  -  r0)  exp(j27ri/0f) 

+  ^2  x\s{t  -  n)  exp (j 2-KUit)  +  n(t),  (2) 

i 

where  bl,To  and  vq  are  the  complex  reflectivity,  delay  and 
Doppler  shift,  respectively,  of  the  target  (if  present),  x\ ,  r, 
and  Vi  are  the  complex  reflectivity,  delay  and  Doppler  shift 
of  the  zth  clutter  scatterer,  respectively,  and  n(t)  is  additive 
noise.  In  (2),  we  have  assumed  that  the  delay  and  Doppler 
shift  of  each  scatterer  remains  constant  across  the  duration 
of  the  Sub-dwell,  which  is  reasonable  if  the  PRF  is  high 
(in  the  order  of  10  kHz).  Note  that  we  make  no  assumption 
on  the  values  of  the  target  or  clutter  Doppler  shift  so  that 
discrimination  of  the  target  based  on  Doppler  differences  alone 
is  not  possible.  We  will  henceforth  assume  a  high  clutter-to- 
noise  ratio  so  that  the  clutter  is  the  dominant  component  and 
additive  noise  is  negligible.  Since  we  only  consider  transmitted 
signals  of  very  short  duration,  the  Doppler  resolution  is  very 
poor.  Therefore,  we  completely  ignore  Doppler  processing  and 
restrict  our  attention  to  delay  or  range  estimation  alone. 

The  received  signal  in  (2)  is  sampled  at  a  rate  fs  to  yield 
the  sequence  gl[n]  =  g\n/ fs),  which  upon  matched  filtering 
yields  the  sequence  yl[n}.  We  define 

as  the  vector  of  matched-filter  outputs  at  delay  or  range  bin 
j,  so  that 


Fig.  1.  Pulse  diagram  of  the  proposed  adaptation  algorithm.  The  pulse 
amplitudes  are  for  illustrative  purposes  only. 


B.  Clutter  and  Target  Models 

Early  investigations  of  waveform  design  for  sea  clutter 
rejection  typically  assumed  a  Gaussian  model  for  the  echoes 
[12],  As  the  spatial  resolution  of  radars  improved,  however, 
this  model  could  not  account  for  the  increased  occurrence 
of  higher  amplitudes  or  spikes.  As  a  result,  a  compound- 
Gaussian  model  for  sea  clutter  has  gained  wide  acceptance 
[14]  and  has  received  theoretical  and  empirical  support.  This 
model  postulates  that  sea  clutter  returns  are  the  result  of 
two  components:  a  speckle-like  return  that  arises  due  to  a 
large  number  of  independent  scattering  centers  reflecting  the 
incident  beam,  and  a  texture  caused  by  large-scale  swell 
structures  that  modulates  the  local  mean  power  of  the  speckle 
return.  Given  the  texture  and  speckle  covariance,  the  complex 
reflectivity  of  the  ith  clutter  scatterer  is  a  circular  complex 
zero-mean  Gaussian  random  vector  [  15]  or  X/  ~  CAf(0, 7j£), 
where  7]  >  0  is  the  texture,  and  £  £  CLxL  is  the  speckle 
covariance  matrix.  The  speckle  remains  correlated  over  a 
short  interval  (~  1 0  ms)  while  the  texture  exhibits  longer 
decorrelation  time  ( ~ 5 0  s)  [17].  Thus  at  the  PRF  we  consider, 
£  ^  II  is  non- white,  where  II/,  is  the  LxL  identity  matrix.  We 
will  assume  that  the  texture  remains  constant  during  a  dwell, 
i.e.  over  a  duration  of  2-3  ms.  We  also  assume  a  Swerling  I 
target  [18],  so  that  b  ~  CA/"(0,  a2 II).  From  (3)  it  follows  that 

y j  ~  CAf(0,a2IK\zs[j  -  n0] |2  +  £/?y). 


N- 1 


yj  =bzs\j  -n0\+  ^2  Zs[n]xj+n,  (3) 

n=-(/V-l) 


where 


b  =  [b0,b\...,bL~']T, 

Xi  =  [x^x],...^^1}71, 

no  is  the  range  bin  that  contains  the  target,  and  zs[m]  = 
s[n]s*[n  —  to],  \m\  <  N  is  the  autocorrelation  function 
of  s[n]  =  s(n/fs),  n  =  0, . . . ,  N  -  1. 


where  the  scalar 

N-l 

(3j  =  ^2  7i+"l2:s[n]l 

n=—(N  —  l) 

is  a  function  of  the  waveform  due  to  its  dependence  on  zs[n}. 

C.  Observations  Model 

At  the  end  of  Sub-dwell  2  of  the  fcth  dwell,  the  measurement 
provided  to  the  tracker  is 

k  [y  fiQ—nv  *  ■  5  y  no )  •  •  •  >  y  ho+nv  ]  > 


where  ho  is  the  range  bin  corresponding  to  the  predicted  target 
range,  nv  =  (14  — 1)/2,  where  14  is  the  number  of  range  bins 
in  the  validation  gate  at  time  step  k  [16],  while  yT  denotes 
the  transpose  of  y.  Since  the  observations  of  target  range  are 
quantized  to  range  bins,  we  will  assume  that  the  actual  target 
range  is  uniformly  distributed  across  the  length  of  the  range 
bin. 

III.  Dynamic  Waveform  Design 

From  (3),  we  can  see  that  the  matched-filter  output  is 
a  convolution  of  the  radar  scene  with  the  autocorrelation 
function  of  the  transmitted  signal.  This  results  in  a  waveform- 
dependent  smearing  of  energy  from  one  range  bin  to  another. 
Clearly,  the  impact  of  the  ‘out-of-bin’  clutter  is  to  reduce  the 
SCR  in  a  particular  range  bin.  Our  objective  is  to  improve 
the  SCR  in  Sub-dwell  2  in  range  bin  ho,  where  the  target  has 
been  predicted.  To  this  end,  we  seek  to  design  S2[n]  so  that 
its  autocorrelation  function  zS2  [to]  is  small  where  the  clutter 
is  strong,  thus  minimizing  the  effect  of  out-of-bin  clutter  in 
the  predicted  target  location.  To  achieve  this,  we  first  require 
estimates  of  the  clutter  power  or  texture  in  the  bins  whose 
scatterers  contribute  to  the  return  y „0. 

A.  Estimation  of  Clutter  Statistics 

Although,  the  texture  is  a  random  process,  the  assumption 
that  it  is  practically  constant  across  a  dwell  permits  us  to  treat 
its  realizations  in  each  range  bin  as  deterministic  but  unknown 
variates.  Accordingly,  suppose  that  we  wish  to  estimate 

®  —  2(N— 1)  >  •  •  •  i  Tn o  5  •  •  •  5  /^~no+2(N—  1)  ? 

From  (3),  we  see  that  there  is  a  many-to-one  linear  mapping 
between  the  scatterer  reflectivities  and  the  matched-filtered 
output  vectors  which  precludes  the  exact  determination  of  0. 
Therefore,  we  instead  seek  its  maximum-likelihood  estimate 
(MLE)  that  maximizes  p( y;  0),  where 

r  T  T  T  l  T 

y  [y no  —  {N—l)'>  *  *  *  Ano’  '  '  •  ^y ho+(N  —  1)J  • 

To  find  this  MLE  however,  we  must  perform  a  complicated 
multi-dimensional  search  since  the  elements  of  y  are  strongly 
correlated.  Instead,  we  use  the  expectation-maximization  (EM) 
[19]  algorithm  to  find  0  that  maximizes  p(x:  0),  where 

_  r  T  T  T  iT 

X  —  [Xn0-2(Ar-l)>  •  •  •  >Xft0’  ’  '  ‘  >  Xfto+2(N-l)J 

is  the  unobserved  or  complete  data  vector.  Starting  with  an 
initial  guess  0^,  the  EM  algorithm  iterates 

U(Q,Gm)  =  £{lnp(x;  ©)|  y,  ©<“>},  (4) 

0(lt+1)  =  argmax(7(0,0(,l)),  (5) 

until  the  change  in  0  falls  below  a  set  threshold.  Note 
that  x  is  complex  Gaussian  distributed  and  its  elements  are 
independent,  given  the  texture  and  speckle  covariance  matrix. 
Together  with  the  fact  that  y  and  x  are  related  through  a  linear 
transformation,  this  simplifies  the  computations  in  (4)  and  (5). 


B.  Phase  modulated  waveform  design 

After  estimation  of  0,  as  described  in  Section  III-A,  let 
us  assume  that  we  have  formed  the  (possibly  disconnected) 
set  ZT  of  range  bin  indices  in  the  Ar  -  1  neighborhood  of 
ho,  where  the  clutter  power  is  large.  For  example,  this  could 
consist  of  the  Nt  <  2(N  —  1)  largest  values  of  the  estimated 
texture  [7^0_(JV_i), . . . ,  tho+{N_1)]  excluding  tho  itself,  as 
it  corresponds  to  the  predicted  target  location.  In  order  to 
minimize  the  out-of-bin  clutter  contribution  to  y n0  in  Sub¬ 
dwell  2,  we  seek  to  minimize  |2S2[to]|2. 

To  achieve  this,  we  choose  a  unimodular  phase-modulated 
waveform  [20], 

s2[n]  =  exp(jA„),  n  =  0, 1, . . . ,  N  -  1, 

as  a  template,  and  seek  to  minimize  the  cost  function 

J(A)  =  XI  1^2  Ml2)  (6) 

where  A  =  [Ao,  Ai, . . . ,  Aat_i]t.  Note  that  the  least  squares 
minimization  in  (6)  is  a  simplified  form  of  a  method  that  may 
be  used  to  design  a  waveform  with  an  arbitrary  autocorrelation 
function.  Further,  if  the  minimization  is  carried  out  over  the 
time  and  frequency  plane,  the  same  approach  can  be  used  to 
design  a  waveform  with  a  specified  ambiguity  function. 

It  is  relatively  straightforward  to  obtain  zS2[m\  [3],  The 
minimization  of  J( A)  in  (6)  can  then  be  easily  accomplished 
using  the  Newton-Raphson  method. 

IV.  Target  Tracking 

The  target  tracker  estimates  the  posterior  probability  density 
function  p(Xfc|Yi;fc)  of  the  target  state  at  dwell  k  given  the 
observations  upto  and  including  dwell  k  [16].  Although  the 
observation  Y provides  a  measurement  of  the  range  of  the 
target,  note  that  it  is  accompanied  by  noise  that  is  uniformly 
distributed  across  the  extent  of  the  range  bin.  Therefore, 
despite  the  fact  that  the  observations  model  is  a  linear  function 
of  the  target  state,  we  cannot  use  the  Kalman  filter  as  the  target 
tracker.  Thus,  we  use  the  particle  filter  [21]  which  estimates 
the  posterior  density  as  a  collection  of  Np  samples  or  particles 
X^.,  i  =1,2,...,  Np,  and  their  corresponding  weights  w\.  At 
each  dwell,  the  density  is  predicted  forward  in  accordance 
with  the  target  dynamics  model  and  updated  according  to  the 
likelihood  function. 

When  targets  are  to  be  tracked  in  the  presence  of  clutter, 
probabilistic  data  association  is  often  used  to  counter  the 
uncertainty  in  the  origin  of  the  measurements  [16].  In  this 
approach,  the  likelihood  function  is  usually  computed  as  the 
average 

p(Yfc|x£)  =  5>(Yfc|x£,nm)p(nm), 

m 

where  Clm  is  an  association  hypothesis.  In  our  case,  however, 
due  to  the  availability  of  estimates  of  the  clutter  and  target 
statistics  in  each  range  bin,  we  are  able  to  directly  compute 
the  likelihood  functions  p( Y^|XJ.). 


Fig.  2.  Magnitude  autocorrelation  function  for  a  linear  frequency  modulated 
(LFM)  chirp  and  the  designed  phase-modulated  waveform.  The  range  bins 
with  strong  clutter  are  marked  with  asterisks. 


Fig.  3.  Comparison  of  averaged  position  RMSE  at  the  end  of  Sub-dwell 
1  (dotted  lines)  and  Sub-dwell  2  (solid  lines)  using  the  designed  waveform. 
The  numbers  on  the  curves  indicate  SCR  in  dB. 


For  this  purpose,  let  nl0  denote  the  range  bin  that  contains 
the  range  corresponding  to  X(. .  With 


=  [zS2  [-(TV  -  1)], . . . ,  *82 [0], . . . ,  zS2  [N  -  1]]T, 


let 


H 


0  zS2  0 


0 


0  Zs 


where  ®  indicates  the  Kronecker  product.  Furthermore,  let 


T  =  diag[7^0  —  nk  j  •  •  •  )  i  ■  ■  ■  i  ^tio+nfe]  ®  53, 
V  —  diag[Vft0  —nv  )  *  •  •  j  VfiQ  ,  .  .  .  ,  Vfi0-\-nv  ]  &  j 


where  n k  =  nv  +  (N  —  1),  Vj  =  &]S[j  —  tiq],  />[•]  is  the 

Kronecker  delta,  and  o2  is  the  MLE  of  a2,  given  that  the 
target  is  in  range  bin  j.  The  covariance  matrix  of  Y /,.,  given 
X.\,  is  then  given  by 

53j  =  H(T  +  V). 


From  (3),  Y&  is  a  zero-mean,  complex  Gaussian  vector  and 
hence  the  weight  w\  oc  p(Yfc|X],)  =  CAf(0,  X,;)  can  be  easily 
computed. 


V.  Simulations 

In  our  simulation  study,  a  target  is  tracked  as  it  moves  in 
the  presence  of  simulated  heavy  sea  clutter.  At  the  start  of 
each  tracking  sequence  of  25  dwells,  the  target  is  located  at 
a  distance  of  10  km  from  the  sensor  and  moves  away  from 
it  at  a  near  constant  velocity  of  5  m/s.  The  clutter  scatterers 
are  distributed  uniformly  in  range  across  the  range  cells,  and 
uniformly  in  Doppler  over  [-1,000,  1,000]  Hz.  The  clutter 

reflectivity  is  sampled  from  a  K-distribution  with  the  necessary 
speckle  and  texture  temporal  correlations  derived  from  real 
sea  clutter  data  as  in  [3],  [17].  The  particles  at  the  start  of 


each  tracking  sequence  are  drawn  fromp(Xo)  =  C7V(Xo,  Po), 
where  Po  =  diag[l,000  10]. 

The  waveform  si[n]  was  chosen  to  be  a  linear  frequency 
modulated  (LFM)  chirp  of  duration  1.5  fis  and  a  frequency 
sweep  of  100  MHz.  We  used  L  =  10  pulses  in  each  sub¬ 
dwell,  a  pulse  repetition  interval  (PRI)  of  100  /is,  and  a  tracker 
update  interval  of  AT  =  1  s.  The  amplitude  of  the  target 
return  was  sampled  from  a  zero-mean,  complex  Gaussian 
process  with  covariance  matrix  ct2Il,  where  a2  was  chosen 
to  satisfy  specified  values  of  SCR.  We  define  the  SCR  to 
be  the  ratio  of  the  target  signal  power  to  the  total  power 
of  the  clutter  in  the  range  bin  containing  the  target.  As  a 
measure  of  the  algorithm’s  performance,  we  present  averaged 
position  RMSE  plots  conditioned  on  convergence,  where  a 
run  was  said  to  have  converged  if  the  true  target  was  always 
located  within  the  validation  gate  [16].  The  plots  are  compared 
with  the  corresponding  performance  obtained  when  the  fixed 
LFM  chirp  waveform  alone  is  transmitted  in  each  dwell,  thus 
representing  a  non-adaptive  system. 

In  Fig.  2,  we  compare  the  magnitude  of  a  typical  phase- 
modulated  waveform,  designed  as  described  in  Section  III 
with  that  of  the  LFM  chirp  transmitted  in  Sub-dwell  1. 
The  range  bins  in  which  the  clutter  was  estimated  to  be 
strong,  or  the  elements  of  the  set  ZT  in  (6),  are  marked  by 
asterisks.  It  can  be  observed  that  in  these  range  bins,  the 
magnitude  autocorrelation  function  is  upto  30  dB  below  the 
corresponding  level  for  the  LFM  chirp. 

Fig.  3  shows  a  comparison  of  the  averaged  position  RMSE 
obtained  at  the  end  of  Sub-dwells  1  and  2,  corresponding  to 
the  transmission  of  the  fixed  LFM  chirp,  and  the  dynamically 
designed  waveform.  We  observe  that  the  RMSE  is  lower  when 
the  dynamically  designed  waveform  is  used  thus  indicating  a 
performance  improvement  of  more  than  20  dB  SCR  for  the 
same  RMSE. 


VI.  Conclusions 

In  this  paper  we  presented  a  waveform-adaptive  tracking  al¬ 
gorithm  that  reduces  tracking  error  by  improving  the  detection 
performance.  Using  a  two-step  approach,  we  first  estimate  the 
clutter  statistics  within  the  validation  gate  and  then  use  these 
estimates  to  design  a  phase-modulated  waveform  for  the  next 
transmission.  The  second  waveform  helps  to  improve  the  SCR 
in  the  target  location  and  thus  improve  detection  performance. 
Since  a  reduction  in  the  number  of  false  alarms  directly 
reduces  uncertainty  in  the  origin  of  the  measurement,  the 
tracking  performance  improves  considerably,  as  demonstrated 
in  a  simulation  study.  Note  that  tracking  performance  can  be 
further  improved  by  adapting  the  waveform  to  the  dynamic 
demands  of  the  tracker  in  addition  to  the  estimates  of  the 
environment.  This  would  permit  tradeoffs  in  the  delay  and 
Doppler  resolution  properties  of  each  waveform  in  addition  to 
their  ability  to  mitigate  sea  clutter. 
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